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ABSTRACT 

Aims. Emission from the 6.7 GHz methanol maser transition is very strong, is relatively stable, has small internal motions, and is 
observed toward numerous massive star-forming regions in the Galaxy. Our goal is to perform high-precision astrometry using this 
maser transition to obtain accurate distances to their host regions. 

Methods. Eight strong masers were observed during five epochs of VLBI observations with the European VLBI Network between 
2006 June, and 2008 March. 

Results. We report trigonometric parallaxes for five star-forming regions, with accuracies as good as ~ 22 yuas. Distances to these 
sources are 2.57+^ kpc for ON 1, 0.776!° ^ kpc for L 1206, kpc for L 1287, 2.38!° ^ kpc for NGC 281-W, and 1.59^ 

kpc for S 255. The distances and proper motions yield the full space motions of the star-forming regions hosting the masers, and we 
find that these regions lag circular rotation on average by ~17 km s _I , a value comparable to those found recently by similar studies. 

Key words. Techniques: interferometric - Astrometry - Masers - Stars: formation - ISM: molecules - Galaxy: kinematics and 
dynamics 



1. Introduction 

Accurate distances and proper motions are crucial for studies of 
the structure and kinematics of the Milky Way. Since we are in 
the Galactic plane, it i s not easy to dete rmine the spiral struc- 
ture of our Milky Way dReid et alj20 09b). Massive star-forming 
regions trace the spiral arms and are objects well-suited to re- 
vealing the structure of the Milky Way. 

Determining the fundamental physical properties of indi- 
vidual objects, such as size scales, masses, luminosities, and 
ages, also depends critically on distance. For example, the dis- 
tance to the Orion Nebula, determined by a number of trigono- 
metric parallax measurements using radi o continuum emis- 
sion from stars in its associated cluster ( Menten e t all |2007t 
Sandstrom et al. 20 07b - water masers dHirota et alj l2007l) . and 
SiO masers dKim et all (2008), turned out to be 10% less than 
previously assumed, resulting in 10% lower masses, 20% fainter 
luminosities, and 20-30% younger ages for the stars in the clus- 
ter. 

The most fundamental and unbiased method of measuring 
distances is the trigonometric parallax, which depends only on 
geometry and is therefore free of any astrophysical assumptions. 
To achieve the sub-milliarcsecond astrometric accuracy at op- 
tical wavelengths requires space-borne observations. The first 
dedicated optical satellite for this purpose was ESA's Hipparcos 
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dPerrvman et alJ[T995l) . With its accuracy of 0.8-2 mas, it was 
capable of measuring distances up to ~200pc - only a small 
fraction of the Milky Way. A new optical astrometry satellite 
GAIA, to be launched in 2012, will be two orders of magni- 
tude more accur ate and push the optical parallaxes to Galactic 
size scales (e.g., iLinde grenl 12009 ) . Despite the high accuracy, 
GAIA will suffer from dust extinction in the Galactic plane, 
in particular in the spiral arms and toward the Galactic cen- 
ter. Here radio astronomy can provide an important comple- 
ment to GAIA, because radio wavelengths do not suffer from 
dust extinction, and Very-Long-Baseline Interferometry (VLBI) 
phase-referencing techniques can provide accuracies better than 
10 yuas, allowing accurate dist ances with errors less than 10% out 
to 10 kpc (see for exam ple iHonma et al]|2007l lHachisuka et al.l 
120091 iReid et"aDl2009ah . 

Strong and compact radio sources, such as molecular 
masers, are ideal targets for trigonometric parallax measure- 
ments. Masers are frequently found in (dusty) star formation re- 
gions (SFRs) and asymptotic giant branch stars. The most nu- 
merous interstellar masers are the 22 GHz water masers, fol- 
lowed in number by the 6.7 GHz methanol masers; to date sev- 
eral h undreds of these masers have been found across the Galaxy 
(e.g.. lPestalozzi et alj|2005t iGreen et al.ll2009t) . They are exclu- 
sively associated with early stages of massive star formation, 
observed both prior to the development of an ultra compact 
Hn (UCHn) region and coexistent wi th these (iMenteni : 1991; 
Ellingsen 2006; Pestalo zzi et al.l f2007). Maser modeling indi- 
cates that the 6.7 GHz emission is likely to be excited via radia- 
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Table 1. Observation information 



Program 


Date 


Participating Antennas^ 


EB32A 


10 Jun 2006 


Jb Wb Ef Mc Tr Nt Hh On Da 


EB32B 


18 Mar 2007 


Jb Wb Ef Mc Tr Nt Hh On Cm 


EB32C 


17 Jun 2007 


Jb Wb Ef Mc Tr Nt Hh On Cm 


GB63A 


31 Oct 2007 


Jb Wb Ef Mc Tr Nt Hh On Yl 


GB63B 


16 Mar 2008 


Jb Wb Ef Mc Tr Nt Hh On Yl 



" The full names of the antennas are: Jb - Jodrell Bank; Wb - 
Westerbork (in single dish mode); Ef - Effelsberg; Mc - Medicina; Tr 
- Torun; Nt - Noto; Hh - Hartebeesthoek; On - Onsala; Da - Darnhall; 
Cm - Cambridge; Yl - single EVLA antenna. 



tive p umping by warm d ust heated by newly formed high-mass 
stars dSobolev et al.ll2007l) . 

Sources of meth anol masers have long lifetimes of ~ 10 4 year 
(Ivan der Waltl 12005). Their velocity spread is typically within 
±5 km s _I about the systematic velocity of the molecular core, 
and single maser features are usually narrow (~lkm s _1 ). The 
kinematics suggest that methanol masers originate from different 
regions than water masers, which have wid er velocity spr eads 
and are produced in protostellar outflows dMentenlll996l) . To 
summarize, the strength, ubiquity, long lifetimes, and small in- 
ternal motions make the 6.7 GHz methanol maser very suitable 
for astrometric purposes. 

Here we present results from observations with the European 
VLBI Network (EVN) of eight strong methanol masers belong- 
ing to well-known massive SFRs in the outer part of the Milky 
Way: ON 1, L 1287, L 1206, NGC281-W, MonR 2, S 252, S 255 
and S269. Recently, H2O maser or 12.2 GHz methanol maser 
parallax measureme nts have been rep orted fo r three of these re- 
gions (NGC 281-W ISato et al.1 120081: S252. iReid etail 120091: 
S 269, lHonma et alj|2007l) . allowing a valuable cross check with 
our measurements. The results of this present work are crucial 
for star formation studies in these regions, and, together with a 
larger sample of parallaxes, will help to understand the structure 
of the Local and Perseus arms. 

In this paper we report on the first parallax measurements 
made with the EVN. Results presented here repla ce prelimi- 
nary, less accurate, and less complete measurements (Ryg l et all 
2008). In the next section we describe the observations and data 
analysis, and in Sect. 3 we explain the method and fitting of the 
parallax and proper motion. The results are then presented per 
source in Sect. 4 and discussed in Sect. 5. The conclusions are 
summarized in Sect. 6. 



2. Observations and data analysis 

The observations were performed with the EVN at five epochs 
between 2006 June, and 2008 March. The exact dates of the ob- 
servations and the participating antennas are listed in Table Q] 
Each observation lasted 24 hours and made use of geodetic-like 
observations t o calibrate the troposphe r ic zenith delays at each 
antenna (see IReid & Brunfhaleri l2004t iBrunthaler et all 120051; 
IReid et al.ll2009al for a detailed discussion). A typical observ- 
ing run consisted of four 6-hour blocks containing ~1 hour of 
geodetic-like observations, ~ 10 minutes of observation of fringe 
finders. The remaining time was spent on maser/background 
source phase-referencing observations. During each run, the 
average on-source time per maser was between ~0.9 and 
~ 1.2 hours depending on the sky position, because half of our 
targets were low declination sources that had a limited visibility. 



Using the technique of phase referencing, each maser was 
observed in a cycle with two (or one) nearby (~1°— 2° separa- 
tion) compact radio quasars, which were used as background 
sources. The sources were switched every 2 minutes. Before 
our EVN observations, we had used the NRAO Very Large 
Array under project AB 1207 in A-configur ation to observe sev- 
eral compact NVSS dCondon et alj IT998) sources within 1° 
from the maser source at two frequencies (8.4 and 4.8 GHz) to 
check their spectral index and compactness. For the best can- 
didates, we determined their position to a sub-arcsecond ac- 
curacy to use them as a position reference in the EVN ob- 
servations. Additionally, several known VLB A c alibrators were 
used namely: J2003+3034 and J0613+1708 (FomalontetaD 



120031) . J0603 +215S(9) and J06 13+1306 dMaetaU 1 19981) 



J0047 +5657 dBeaslev et alJ l2002h J2223+6249 dPetrov et alJ 
120051). J0035+ 6130 dPetrovetalJ [2006b . and J0606-0724 
dKovalev et alj[2007b . Table [2] lists the positions of the masers 
and their background sources. 

The observations were performed with eight intermediate 
frequency bands (IF) of 8 MHz width, each in dual circular po- 
larization sampled at the Nyquist rate with 2 bits per sample, 
yielding a recording rate of 512 Mbps. The data were correlated 
in two passes at the Joint Institute for VLBI in Europe (JIVE), 
using an integration time of 0.5 seconds, affording a field of 
view of 1/2 (limited by time-averaging smearing). The maser 
data were correlated using one 8 MHz IF band with 1024 spec- 
tral channels, resulting in a channel separation of 7.81kHz or 
0.41 km s" 1 at 6.7 GHz. The quasar sources were correlated in 
continuum mode with eight IFs of 8 MHz width with a channel 
separation of 0.25 MHz. 

The data were reduced using the NRAO Astronomical Image 
Processing System (AIPS). The geodetic-like observations were 
reduced separately, and tropospheric delays were estimated for 
each antenna. The data were reduced following the EVN guide- 
lines, applying parallactic angle and ionospheric delay correc- 
tions. The JIVE correlator model uses Earth's orientation pa- 
rameters, which are interpolated from the appropriate daily- 
tabulated values, so it is not necessary to correct them after the 
correlation. The ionospheric delays were based on the JPL GPS- 
IONEX total vertical electron content maps of the atmosphere. 
Amplitudes were calibrated using system-temperature measure- 
ments and standard gain curves. A "manual phase-calibration" 
was performed to remove delay and phase differences between 
the IFs. The Earth rotation was corrected for with the task 
'CVEL. For each maser, a spectral channel with one bright and 
compact maser spot was used as the phase reference. The data 
was Hanning-smoothed to minmize Gibbs ringing in the spectral 
line data. To avoid the strong fluctuations caused by the band- 
pass edges, the outer two ch annels in each IF were discarded 
(following IReid et alJ.l20 09a). The positions of the masers and 
background sources were extracted by fitting 2D Gaussians to 
the maps. 

Here, we report the results on the first five sources for 
which we have completed the analysis: ON 1, L 1206, L 1287, 
NGC 281-W, and S 255. For the other sources, S 252, S 269, and 
MonR 2, we had problems with the calibration of the data, which 
were likely caused by residual atmospheric delay, so that we 
were unable to achieve sufficient accuracy for a parallax mea- 
surement. 

3. Method and fitting 

The absolute positions of the masers were not known with milli- 
arcsecond accuracy before our first observation. Only after the 
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Source 


R.A. (J2000) 


Dec. (J2000) 




P.A." 


Brightness^ 


Restoring 


Beam* 




(h:m:s) 


(° : ' : ") 


(°) 


C) 


(Jy beam" 1 ) 


(mas, mas 


, deg) 


ON 1 


20:10:09.074 


+31:31:35.946 






1-7 


5.7 x 3.9, 


86 


J2003+3034 ... 


20:03:30.244 


+30:34:30.789 


1.71 


-124 


0.124 


5.3x3.7, 


76 


J2009+3049 ... 


20:09:17.588 


+30:49:24.580 


0.73 


-14 


0.008 


5.6x3.7, 


88 


L 1206 


22:28:51.407 


+64:13:41.314 






0.3-1.4 


5.7 x 4.0,- 


-80 


J2223+6249 ... 


22:23:18.097 


+62:49:33.805 


1.53 


-157 


0.093 


5.7x3.9,- 


-84 


J2225+6411 ... 


22:25:27.993 


+64:11:15.030 


0.37 


-96 


0.002 


5.6x3.8,- 


-80 


L 1287 


00:36:47.353 


+63:29:02.162 






0.5-4.1 


6.0x4.1,- 


-57 


J0035+6130 ... 


00:35:25.311 


+61:30:30.761 


1.98 


-176 


0.092 


6.0x4.1,- 


-58 


J0037+6236 ... 


00:37:04.332 


+62:36:33.310 


0.88 


178 


0.061 


5.9 x 4.2,- 


-57 


NGC281-W ... 


00:52:24.196 


+56:33:43.175 






1-12 


5.7 x 4.2,- 


-56 


J0047+5657 ... 


00:47:00.429 


+56:57:42.395 


0.84 


62 


0.142 


5.6 x 4.0,- 


-65 


J0052+5703 ... 


00:52:54.303 


+57:03:31.460 


0.50 


-8 


0.010 


5.6 x 3.9,- 


-62 


MonR2 


06:07:47.862 


-06:22:56.518 






9-22 


8.6 x 4.4, 


42 


T0606-0724 


06-06-43 546 


-07-24-30 232 


1.06 


-165 


0.135 


9.0 x 4.1, 


40 


S 252 


06:08:53.344 


+21:38:29.158 






17-44 


6.4x3.9, 


56 


J0603+215S .. 


06:03:51.557 


+21:59:37.698 


1.22 


-73 


0.019 


6.7x3.8, 


54 


S 255 


06:12:54.020 


+ 17:59:23.316 






4 


6.9 x 3.9, 


56 


J0613+1708 ... 


06:13:36.360 


+ 17:08:24.946 


0.87 


169 


0.040 


7.1 x3.9, 


46 


S 269 


06:14:37.055 


+13:49:36.156 






0.2-1.1 


6.8x3.9, 


58 


J0613+1306 ... 


06:13:57.693 


+ 13:06:45.401 


0.73 


-167 


0.087 


6.8x3.8, 


57 



" Separation, (f>, and the Position Angle (east of north), RA. , between the maser and the background source. 
h The brightness and restoring beam (east of north) are listed for the second epoch. 



first epoch was analyzed were accurate positions determined. As 
a result, different correlator positions were used between the first 
and following epochs. In the first epoch, position errors above 
0."2 made it difficult to calibrate the phases because of a high 
fringe rate, which could not be interpolated correctly when ap- 
plied to the data. Shifting the maser into the phase center turned 
out to be very time-consuming with available software because 
AIPS and the EVN correlator use apparent positions calculated 
at different times (AIPS at UT of the first observing day versus 
the EVN correlator at the end of the first scan). Even an infinites- 
imally small shift in the epoch coordinate using the AIPS tasks 
'CLCOR' or 'UVFIX' can result in a substantial position change 
and phase shifts from the old EVN apparent position to the new 
AIPS apparent position at UT. In a parallax experiment, accu- 
rate and consistent positions throughout the experiment are cru- 
cial. Since it was difficult to calibrate the phases for masers with 
a large position offset or to shift these to the phase center, several 
sources in the first epoch could not be used for phase referenc- 
ing. Although it should be possible to correct for this difference 
in registration, it would have required a software effort beyond 
the scope of the current project. Therefore we chose to discard 
some of the first-epoch data. 

Most of the ionospheric delay is removed using measure- 
ments of the total electron content of the ionosphere obtained 
from dual-band global positioning system measurements (e.g., 
iRos et al1l2000l) . The expected residual delay is still a few cen- 
timeters, comparable to the residual tropospheric delay at the 
observing frequency. Our single frequency geodetic-like obser- 
vations do not allow a separation of the ionospheric and tro- 
pospheric contributions. The measured group-delays from the 
geodetic blocks are interpreted as having a tropospheric origin, 
and the data are phase-delay-corrected accordingly. An iono- 
spheric group delay makes a different contribution to the phase 
delay, so that the tropospheric correction can even deteriorate 



the calibration when the group delay has a partially ionospheric 
origin. Since the residual tropospheric delay is a few centimeters 
at 6 GHz and can easily be confused with an ionospheric resid- 
ual, the tropospheric correction is not expected to be very useful 
in improving the quality of calibration. As a test, all data were 
reduced with and without a "tropospheric" delay correction. In 
slightly more than half of the cases, the signal-to-noise ratio im- 
proved by applying the tropospheric correction. We noted a trend 
toward improvement with decreasing declination: both high dec- 
lination sources, 5 - 64°, L 1206 and L 1287 improved in 2 
out of 5 epochs; NGC281-W at 5 57° improved in 3 out 
5 epochs; and ON 1 at 5 - 31° even improved in 4 out of 5 
epochs. Ionospheric delay saturates at zenith angles greater than 
~70°. However, the tropospheric delay continues to grow rapidly 
at larger zenith angles and can dominate the ionos pheric delay at 
zenith angles of ~ 80° (see Thompson et al. 199 if. 

For each SFR, we found the emission to arise from a num- 
ber of separate maser spots. Most of the line profiles stretched 
over several channels. Both spatially and/or spectrally different 
maser spots were considered as distinct maser components. We 
inspected the behavior of the proper motion for each maser spot 
relative to the reference spot. Maser spots with strong nonlinear 
proper motions or a large scatter of position about a linear fit 
were discarded. Only compact maser spots with well-behaved 
proper motions were used for the parallax fitting. 

The average internal proper motion of the maser spots ranged 
between 0.06 and 0.23 mas yr . Considering the distance of 
each maser, these proper motions correspond to 0.5-1 km s , 
much lower than the internal proper m otions of water masers , 
which can reach up to 20-200 km s -1 (Hachis uka et al.l 12006). 
The only exception was ON 1, which separates into two distinct 
maser groups with a relative proper motion of 0.52 mas yr -1 , or 
6.3 km s . This particularity is discussed in Sect. 3.1. 
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The parallaxes and proper motions were determined from 
the change in the positions of the maser spot(s) relative to the 
background source(s). The data were fitted with a parallax and 
a linear proper motion. Since the formal position errors are only 
based on the signal-to-noise ratios determined from the images, 
they do not include possible systematic errors from residual 
zenith delay errors or source structure changes. This leads to a 
high reduced x 1 value for the fits, so we added error floors in 
quadrature to the positions until reduced^ 2 values close to unity 
were reached for each coordinate. 

First, we performed parallax and proper motion fits for each 
maser spot relative to one background source. Then we made 
combined fits with respect to each background source, assuming 
one parallax but different proper motions for each maser spot. 
Finally, we repeated this combined fit for both the background 
sources together. The position measurements of different maser 
spots are not independent, since systematic errors, such as an 
unmodeled atmospheric delay, will affect all maser spots in a 
similar way. If these systematic errors dominate, this will lead 
to unrealistically small errors. The most conservative approach, 
which we adopted, is to assume that the systematic errors are 
100% correlated. Then the error of the combined fit has to be 
multiplied by yN, where N is the number of maser spots. 

However, this will overestimate the error, if significant ran- 
dom errors are present (e.g., owing to maser blending and struc- 
tural changes over time), since the latter are not correlated be- 
tween different maser spots. Random errors can be reduced by 
averaging the p ositions of the different maser spots (followin g 
the approach of lBartkiewicz et alj|2008t lHachisuka et alj |2009). 
We calculated the average positions with respect to each back- 
ground source after removing their position offsets and proper 
motions. Then, we performed a parallax fit on these averaged 
data sets relative to each individual background source, and on 
both the background sources combined. This approach has the 
advantage that we can reduce the random errors, while leaving 
the systematic errors intact. In Table[5]we list the individual par- 
allax and proper motion fits to each maser spot, the combined fit 
of all the maser spots and the fit of the averaged data sets. 

For most masers, we observed two background sources. 
Some of these background source pairs show a variation in their 
separation (mimicking a proper motion) of up to 1-2 mas yr _1 
(on average this apparent movement was 0.54 mas yr _I ). This 
apparent movement is much greater than expected for extra- 
galacti c sources measured at 12 GHz, < 0.1 mas yr~' dReid et al.l 
2009a) or 22 GHz, < 0.02 mas yr" 1 dBrunfhaler et alj 120071) . 
While at higher frequencies, the radio emission is typically dom- 
inated by the flat spectrum radio cores, the steep spectrum emis- 
sion from the radio jets becomes stronger at lower frequencies. 
Thus, structure changes in these jets can lead to apparent mo- 
tions, even if these jets are unresolved. At the distances of the 
respective masers, these apparent movements correspond to 2— 
10 km s . Since we do not know which of the two background 
sources is responsible for this movement, the apparent move- 
ments of the background source pairs were added in quadrature 
to the final errors of the averaged proper motions of the masers. 
A summary of our parallaxes and proper motions is given in 
Tabled 

4. Individual sources 

4.1. Onsala 1 

The maser emission of Onsala 1 (ON 1) consists of two groups 
separated spatially by ~940 mas, which have different radial ve- 



Table 3. Parallax and proper motion results 



Source 


71 


Ar 




Us 




(mas) 


(kpc) 


(mas yr" 1 ) 


(mas yr -1 ) 



ON 1 0.389 ± 0.045 2.57+°^ -3.24 ± 0.89 r -5.42 ± 0.46^ 

L1206 1.289 ± 0.153 0.776^ 0.27 ± 0.23 -1.40 ± 1.95 
L1287 1.077 ± 0.039 0.929^4 -0.86 ± 0.11 -2.29 ± 0.56 
NGC281-W 0.421 ±0.022 2.38*j;}| -2.69 ± 0.16 -1.77 ± 0.11 
S255 0.628 ±0.027 1.59+°°' -0.14 ± 0.54 -0.84 ± 1.76 



" The proper motion, /j„, includes the factor cos(<5). 

* The errors of the proper motion take the uncertainty resulting from 
an apparent nonzero movement between the background sources into 
account. 

c For the proper motion of ON 1, we took an average of the north and 
south components. 

15 — — 




25 20 15 10 5 -5 -10 



East Offset (mas) 

Fig. 2. Phase-referenced images for two background sources be- 
longing to maser ON 1, J2003+3034 (top) and J2009+3049 (bot- 
tom) in the third epoch. Position offset (0,0) corresponds to the 
position listed in Table [2] The contour levels start at a 3<x level, 
6.1, and 0.6 mJy beam -1 , respectively, and increase by factors of 

VI 



locities by 15kms _1 (Fig. [TJ. The northern group, which in- 
cludes the reference channel, is centered at a vlsr of * km s _1 
and the southern group at a; 15 km s . Four maser spots in the 
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Fig. 1. Velocity-integrated maps of the northern (top, left) and southern (bottom) maser groups in ON 1, together with the spectrum 
(top, right). Position offset (0,0) corresponds to the position listed in Table[2] Maser spots are indicated with color codes for different 
radial velocities. The areas of the colored circles and the colored histogram entries in the spectrum are scaled to the peak flux of that 
spot. The black line in the spectrum is the intensity of the maser within a selected surface, which is not necessarily the same as the 
intensity of the maser spot retrieved from a Gaussian fit. For the northern group, the contour levels start at 0. 1 Jy beam 1 km s , in 
the southern group 1.0 Jy beam -1 km s _1 , and increase by factors of V2. 



northern group were suitable for parallax fitting. The masers in 
the southern group were not used for the parallax fit, since the 
phase-calibration was less accurate because of the large distance 
to the phase-reference center (located in the northern maser 
group). We find a parallax of 0.389 ± 0.045 mas, corresponding 
to a distance of 2.57^2^ kpc. The proper motions of the south- 
ern maser spots were fitted by assuming the parallax result of the 



northern group. All results are listed in Table [5] and the parallax 
fit is plotted in Fig. [3] 

For ON 1, two background sources were detected in epochs 
two through five. The last (fifth) epoch had poor (u, v) cover- 
age, which resulted in very elongated synthesized beams in the 
east-west direction (J2003+3034, 12.1 x 3.7 mas 2 , P.A.~ 87°; 
J2009+3049, 1 1.0x3.3 mas 2 , P.A.~ 88°) compared to the repre- 
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Fig. 3. Results of the parallax fit based on four maser spots. The left graph shows a combined fit on all data, the right graph is a fit 
on the averaged data sets. The filled colored dots mark the data points in right ascension, while the filled colored triangles mark the 
declination. The solid line is the resulting fit in right ascension, the dashed line in declination. Different colors indicate a different 
background source. 
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Fig. 4. Variation of the separation between background sources 
J2003+3034-J2009+3049 belonging to ON 1. The solid line and 
the dots represent the right ascension data, while the dashed line 
and the filled triangles represent the declination data. 



4.2. L 1206 

In L 1206, two maser groups were found separated by ~ 100 mas. 
A third maser group was found northeast of the other two. This 
third maser group was not detected in all the epochs because of 
its weakness so was not used in the parallax measurements. Also 
the maser spot at - 1 1 .6 km s -1 was omitted for the same reason. 
The masers in L 1206 are shown with a spectrum in Fig. [5] The 
parallax fit used four compact maser spots and resulted in 1 .289± 
0.153 mas, corresponding to a distance of 0.776^y ( 'g3kpc. The 
results of the parallax and proper motion fits are displayed in 
Fig.Qand listed in Table[5] 

Only one of the two background sources, J2225+6411, 
was detected in the first epoch. The other background source, 
J2223+6249, had a much larger separation from the maser (i.e., 
the phase reference) and the transfer of the phase solutions to 
the J223+6249 data probably failed. The apparent movements 
between the background sources were 0.22 + 1.45 mas yr~' in 
right ascension and 1.91 + 1.43 mas yr~' in declination (Fig. [8j. 



sentative numbers listed in Table [2] As a result, the right ascen- 
sion data had a large uncertainty in epoch five. We estimated this 
uncertainty from fitting the variation in the background-source 
pair separation to be 2 mas (see Fig.[4]i. Since this was a com- 
bined error for both the background sources, the position error 
for each individual background source was V2 mas. We added 
an additional error floor of V2mas to the error given for all 
the right ascension data points in the fifth epoch. The apparent 
movements between the two background sources, J2003+3034- 
J2009+3049, were 0.88 ± 0.06 mas yr~' in right ascension and 
0.45 + 0.21 mas yr in declination (Fig. [4]). 



4.3. L 1287 

Toward L 1287 we find three maser groups, two close together 
(within ~ 40 mas) and the third ~ 1 35 mas southward, shown 
with a spectrum in Fig. [9] A total of six maser spots from the 
three groups were used for the parallax fit. We report a paral- 
lax value of 1 .077 + 0.039 mas, corresponding to a distance of 
0.929^Q33kpc. The results are plotted in Fig. QT|and listed in 
Table |5] 

The variation of the separation between the two background 
sources, J0035+6130 and J0037+6236, is shown in Fig.[[I] Two 
data points appear to be outliers: the right ascension at epoch 
four and declination at epoch two. Increasing the data errors to 
0.9 and 0.7 mas, respectively, for these outliers yields an appar- 
ent movement fit of -0.05 + 0.12 mas yr -1 in right ascension 
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Fig. 7. Results of the parallax fit for L 1206 based on four maser spots. The left graph shows a combined fit on all data, the right 
graph is a fit on the averaged data sets. The filled dots mark the data points in right ascension, while the filled triangles mark the 
declination. The solid line is the resulting fit in right ascension, the dashed line in declination. Different colors indicate different 
background sources. 
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Fig. 11. Results of the parallax fit for L1278 based on six maser spots. The left graph shows a combined fit on all data, while the 
right graph is a fit on the averaged data sets. The filled dots mark the data points in right ascension, while the filled triangles mark the 
declination. The solid line is the resulting fit in right ascension, the dashed line in declination. Different colors indicate a different 
background source. 



and 0.54 + O.lOmasyr -1 in declination (Fig. [T2l. We investi- 
gated the visibility amplitudes as a function of (u, v) distance 
and found evidence of extended structure. For J0037+6236 at 
the fourth epoch, evidence of extended structure in the right as- 
cension direction was found; however, for J0035+6130 at the 
second epoch, we found no indications of such structure. 



4.4. NGC281-W 

We found four groups of methanol masers toward NGC281- 
W, as illustrated in Fig. [T3l together with the source spectrum. 
For the parallax fit we used a total of six maser spots coming 
from three groups, the central, southern, and southeast groups 
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Fig. 15. Results of the parallax fit for NGC 281-W based on six maser spots. The left graph shows a fit on all data, while the right 
graph is a fit on the averaged data sets. The filled dots mark the data points in right ascension, while the filled triangles mark the 
declination. The solid line is the resulting fit in right ascension, the dashed line in declination. Different colors indicate a different 
background source. 
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Fig. 16. Left: Variation of the separation between background sources J0047+5657-J0052+5703 belonging to NGC 281-W. Right: 
the proper motion fit on the component of J0052+5703 with respect to the central source J0052+5703. The solid line and the dots 
represent the right ascension data, while the dashed line and the filled triangles represent the declination data. 



(Fig. [T3T >. The weaker maser spots, like the ones located in the 
western group, were omitted in the parallax fitting. We find 
a parallax of 0.421 ± 0.022 mas corresponding to distance of 
2.38!° !^ kpc. The results are shown in Fig.[l5]and listed in Table 

For NGC-281-W, two background sources were detected 
successfully at epochs two through five. Background source 
J0052+5703 showed a double structure with the weaker com- 
ponent located ~25mas westward of the stronger component 
(Fig. [T4l >. The proper motion of the western component with 
respect to J0052+5703 was jj a = 0.64 + 0.63 mas yr" 1 and 
jjg = -0.36 + 0.39 mas yr (Fig. [TBI , i.e, exhibiting no de- 
tectable proper motion. The apparent movement between the 
two background sources, J0047+5657 and J0052+5703, was 
also close to zero: -0.16 ± 0.13 mas yr~' in right ascension and 
-0.1 1 ±0.11 mas yr~' in declination. 



4.5. S255 



For S 255, we could only use emission from one maser chan- 
nel, namely the reference channel. Based on this maser spot, at 
4.6 km s _1 (Fig.fTTli. we find a parallax of 0.628 ±0.027 mas cor- 
responding to a distance of 1.59^ ^ kpc. The results are shown 
in Fig.[T9land listed in Table[5] 



We observed only one background source, J0613+1708, for 
S 255 (Fig.fTSTl, which was detected in epochs two to five. With 
only one background source, we could not check that its appar- 
ent movement was small, and so, conservatively, we added un- 
certainties of 1 mas yr~' in both coordinates. 
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Fig. 5. Velocity-integrated map and spectrum of L 1206. 
Position offset (0,0) corresponds to the position listed in Table 
|2] Maser spots are indicated with color codes for different radial 
velocities. The areas of the colored circles and the colored his- 
togram entries in the spectrum are scaled to the peak flux of that 
spot. The black line in the spectrum is the intensity of the maser 
within a selected surface, which is not necessarily the same as 
the intensity of the maser spot retrieved from a Gaussian fit. The 
diffuse and weak spots were omitted, such as the spot at (-10,95). 
Contour levels start at 0. 1 Jy beam -1 km s _1 and increase by fac- 
tors of V2. 
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Fig. 6. Phase-referenced images for two background sources be- 
longing to maser L 1206, J2223+6249 (top) and J2225+6411 
(bottom) in the fifth epoch. Position offset (0,0) corresponds to 
the position listed in Table [2] Contour levels start at a 3<x level, 
3.5, and 0.5 mJy beam -1 respectively, and increase by factors of 
V2. 



5. Discussion 

5.1. Space motions 

With parallax and proper motion measurements, one can calcu- 
late the full space motion of the masers in an SFR in the Galaxy. 
Using an accurate model of Galactic dynamics (Galactic rotation 
speed, ©o, distance of the Sun to the Galactic center, Rq, rota- 
tion curve) for a source and removing the modeled contribution 
of Galactic rotation allows one to retrieve source peculiar mo- 
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Fig. 8. Variation of the separation between background sources 
J2223+6249-J2225+6411 belonging to L1206. The solid line 
and the dots represent the right ascension data, while the dashed 
line and the filled triangles represent the declination data. 



tion relative to a circular orbit. The peculiar motion is described 
by the vectors U, V, and W, locally toward the Galactic center, 
toward the direction of rotation and toward the North Galactic 
Pole, respectively. 

In Table HI the peculiar motions of our sources are given 
for two different Galactic models: the I.A.U. recommended val- 
ues for the LSR motion Rq = 8.5 kpc, ©o = 220 km s _1 in 
columns 2-4; a new model with Rq = 8.4 kpc, ©o = 254 km s , 
based on a large sample of trigonometric parallax measure- 
ments dReid et al.ll2009bl) . in columns 5-7. Both models adopt 
the H ipparcos values for the Solar motion (Dehnen & Binney 
1 19981) and assume a flat rotation curve. We used the parallaxes 
and proper motions as are listed in Table[3] where the uncertainty 
introduced by the apparent movement of the background sources 
was included in the proper motion uncertainty of the source. 

In Fig. [20] we plot the SFRs studied in this work at their 
determined distances with their peculiar motion in the Galactic 
plane, after removing the Galactic rotation. For either model, it 
is clear that three SFRs show a s i milar la g on circular rotation 
as the SFRs studied by lReid et al.1 (l2009bl) . while two have a cir- 
cular rotation close to zero lag. For NGC281-W, the zero lag 
can be understood from its large distance from the Galactic disk 
(b - -6°) and the dominant contribution from the expanding su- 
per bubble to the peculiar motion of this SFR (see Sect. 15. 5b . For 
S 255, it is not clear why its circular velocity is close to Galactic 
rotation (Sect. I5.61 >. The motion towards the Galactic center, U, 
tends to be positive for most of the SFRs. 

Recently, Mc Millan & Binneyl d2009h have publ ished a re- 
analysis of the maser astrometry presented by Rei d et ail 
(2009b). They fitted the data with the r evised solar peculiar 
motion of Binney (2009), mentioned in Mc Millan & B inney 
(2009), where the V component of the solar peculiar veloc- 
ity, V©, is increased from 5.2 km s _1 to 11 km s _1 . With a 
larger V the lag of SFRs on Galactic rotation would decrease 
to ~ 1 1 km s , which agrees with the expected velocity disper - 
sion for young stars of ~ 10 km s _1 dAumer & Binnevl l2009). 
However, even when considering the revised V , we note that 
most SFRs still rotate more slowly than the Galactic rotation 
(by ~ 1 1 km s -1 ) and that their velocities are not randomly dis- 
persed. 
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Fig. 9. Velocity-integrated map and spectrum of L 1287. Position 
offset (0,0) corresponds to the position listed in Table [2] Maser 
spots are indicated with color codes for different radial veloci- 
ties. The areas of the colored circles and the colored histogram 
entries in the spectrum are scaled to the peak flux of that spot. 
The black line in the spectrum is the intensity of the maser within 
a selected surface, which is not necessarily the same as the in- 
tensity of the maser spot retrieved from a Gaussian fit. Contour 
levels start at 0.05 Jy beam -1 km s _1 and increase by factors of 
V2. 



5.2. Onsala 1 

Our trigonometric parallax measurement places ON 1 on a dis- 
tance from the Sun of 2.51^^ kpc, somewhat closer than the 
(near) kinematic distance of 3.0 kpc, based on the methanol 
maser line at 15 km s -1 . The latter assumes a flat rotation 
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Fig. 10. Phase-referenced images for two background sources 
belonging to maser L 1287, J0035+6130 (top) and J0037+6236 
(bottom) in epoch three. Position offset (0,0) corresponds to the 
position listed in Table[2] Contour levels start at a 3cr level, 3.6, 
and 2.4 mJy beam -1 respectively, and increase by factors of y2. 
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Fig. 12. Variation of the separation between background sources 
J0035+6130-J0037+6236 belonging to L1287. The solid line 
and the dots represent the right ascension data, while the dashed 
line and the filled triangles represent the declination data. 



curve with as distance of the Sun to the Galactic center, Rq = 
8.5 kpc, and for the Galactic rotation speed, ©o = 220 km s _1 . 
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Fig. 13. Velocity-integrated map and spectrum of NGC 281-W. 
Position offset (0,0) corresponds to the position listed in Table 
|2] Maser spots are indicated with color codes for different ra- 
dial velocities. The areas of the colored circles and the colored 
histogram entries in the spectrum are scaled to the peak flux 
of that spot. The black line in the spectrum is the intensity of 
the maser within a selected surface, which is not necessarily the 
same as the intensity of the maser spot retrieved from a Gaussian 
fit. Contour levels start at 0.5 Jy beam -1 km s -1 and increase by 
factors of y2. 



However, it is farther than the commonly adopted (near) kine- 
matic distan ce of 1.8 kpc based on formaldehyde at v'lsr = 
11.2km s" 1 dMacLeodetal.H l998) using the Galactic rotation 
curve of W outerloot & Brand! d 19891) . with R = 8.5 kpc, ® = 
220 km s -1 . Our result resolves the near/far kinematic distance 
ambiguity and locates ON 1 in the Local spur consistent with the 
near kinematic distance. 

We found different proper motions for the northern and 
southern maser groups (Table [5j. The resulting difference in 
the east-west direction was small, fx a = 0.18 + 0.24 mas yr -1 , 
but in the north-south direction, the southern group was moving 
by us = -0.77 ± 0.12 mas yr -1 away from the northern group. 
At a distance of 2.6 kpc, this corresponds to a relative speed of 
9.4km s -1 . 

This large proper motion may be explained by the masers 
being located in the molecu lar gas surrounding an e xpanding 
Hn re gion, as suggested by iFish & Reidl d2007l) and ISu et al.l 
(2009). The Hn region is located at {a = 20 h 10 m 09. s 03, 6 = 
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Fig. 14. Phase-referenced images for two background sources 
belonging to maser NGC281-W, J0047+5657 (top) and 
J0052+5703 (bottom) in epoch five. Position offset (0,0) corre- 
sponds to the position listed in Table [2] Contour levels start at 
a 3cr level, 2.8, and 0.6 mJy beam -1 respectively, and increase 
by factors of y2. The bottom image shows also the additional 
component ~ 25 mas westward of J0052+5703. 



+3r31'35."4, J2000) and between the two maser groups. The 
radial velocity of the Hn region, from th e H76a recombination 
line, is 5. 1+2.5 kms" 1 dZheng et all 19851) . The methanol masers 
would be (with respect to the rest frame of ON 1 at 5.1 km s -1 ) 
in a blue-shifted component at ~0 km s _1 , northward, and a red- 
shifted component at ~15 km s -1 , southward of the Hn region. 
Also, the hydro xyl masers in ON 1 provide some confirmatio n 
of this scenario dNammahachak et alj|2006t iFish & Reidll2007l) . 
The expansion velocities of the masers, calculated by assuming 
a linear expansion from the center of the Hn region, would be 
5.8 km s -1 and 3.6 km s -1 for the northern and southern maser 
groups, respectively. 

5.3. L 1206 

For L 1206, the trigonometric parallax distance, 0.116^^* kpc, 
is shorter than the kinematic distance; 1 .0 kpc (based on a hy- 
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Fig. 17. Phase-referenced image for maser S 255 at channel 4.6 
km s _1 . Position offset (0,0) corresponds to the position listed 
in Table |2] Contour levels start at 0.23 Jy beam -1 , and increase 
by factors of ^[2. 
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Fig. 18. Phase-referenced images for background source 
J0613+1708 belonging to maser S255, in the third epoch. 
Position offset (0,0) corresponds to the position listed in Table 
|2] The contour levels start at a 3<x level, 3.0 mJy beam -1 , and 
increase by factors of y2. 



droxyl maser, at vlsr = ~8.5 km s 1 , us ing the Galactic rota- 
tion curve of rWouterloo t & Brand! dl 989ft with R = 8.5 kpc, 
0o = 220 km s~\ iMacLeod etall 1 19981) . 1.23 kpc (ammo- 
nia, vlsr = ~9-9 km s _1 , using the Galactic rotation curve of 
iBrand & Blitz] d 19931) . with R = 8.5 kpc, O = 220 km s -1 
iMolinari etall 1 1 996h and 1.4 kpc (6.7 GHz methanol, flat ro- 
tation curve, Ro = 8.5 kpc, Oo = 220 km s -1 ). Our result 
places L 1206 in the Local spur. Cep A, a nearby (3.°8 sep- 
aration, PA. 126 °) star-forming region, is located at a dis- 
tance of 700 ± 39 pc, as determined by a parallax measure - 
ment of 12.2 GHz methanol masers (Mos cadellietall [2009). 
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Fig. 19. Results of the parallax fit for S 255 based on one maser 
spot at 4.6 km s . The filled dots mark the data points in right 
ascension, while the filled triangles mark the declination. The 
solid line is the resulting fit in right ascension, the dashed line in 
declination. 



The L 1206 SFR therefore seems to be in the same part of 
the Local arm as Cep A. L 1206 is a dark cloud with an in- 
frared source, IRAS 22272+6358A, which coincides with the 
methanol maser em ission and a 2.7 mm dust continuum peak 
(Belu-anetai]|2006). Since neither 2 cm nor 6 cm radio emis- 
sion is detected dWilking et alJll989t iMcCutcheon et al.lll99lT) . 
L 1206 is thought to be in a young phas e prior to the forma- 
tion of an Hn region. Beltra net al.l d2006l) report large CO out- 
flows, and put the systemic velocity of the ambient medium be- 
tween [-13.5, -8.5] km s , consistent with our systemic veloc- 
ity range of [-13.3, -10.9] km s _1 . 

5.4. L 1287 

The parallax for L 1287 sets it at a distance of 0.929+°°|* kpc. 
This is close to the photometric distance of ~850pc from 
dYang et alj 1 19911) . which placed L 1287 in the Local arm. 
However, the kinematic distance based on the methanol maser 
line would place L 1287 at 2kpc in the Perseus arm. 

Methanol masers in the dark cl oud L 1287 are located at 
the base of the bipolar CO outflow (Yang et al. 1991) originat- 
ing in the main core of the cloud. The infrared point source 
(IRAS 00338+6312) in the center of the core indicates a proto- 
stellar object surrou nded by cold, high density NH3 (1,1) gas 
(lEstalella et alJfl993l) . 

5.5. NGC281-W 

There are numerous distance estimates for NGC281-W in 
the literature. Opt ical photometry places the cloud between 
2 and 3.68 kpc flSMmlessl 1 1954 ; ICruz-Gonzalez et all 1 1974 
lHenningetal.1 11994 iGuetter & Turnerlll997l) . while the kine- 
matic distance is 3 kpc (iLee & Juna l2003). bas ed on a vlsr 7 
-30 km s _1 , using the Galactic rotation curve of lClemensI (1985) 



[axes of 6.7 GHz methanol masers 13 




Fig. 20. An artist's impression of a plane-on view of our Galaxy 
(image credit: R. Hurt NASA/JPL-Caltech/SSC), overlaid with 
parallax measurements of wa ter and methanol m asers between 
60° < I < 240° taken from iReid etaT] d2009bl) . in which re- 
cent parallax measurements have been put together to study the 
Galactic structure. The black dot marks the Galactic center, the 
yellow dot encircled with black the Sun, the yellow dots mark 
the parallaxes taken from the literature, the orange dots encir- 
cled with white, and white arrows mark the parallaxes and pe- 
culiar motions obtained in this work. The peculiar motions are 
shown after removing the Galactic rotation assuming the values 
of R = 8 .4 kpc and O = 254 km s" 1 . 



and /?o = 8.5 kpc, 0q = 220 km s -1 . The radial velocity of the 
methanol maser line at -29 km s _1 , assuming a flat Galactic ro- 
tation curve, and Rq = 8.5 kpc, So = 220 km s _1 , would place 
NGC281-W at2.5 kpc. 

Recently, ISato et all d2008l) reported a water maser paral- 
lax value of 0.355 ± 0.030 mas, measured with the VERA in- 
terferometer, corresponding to a distance of 2.82 + 0.24 kpc. 
We measured a parallax of 0.421 ± 0.022 mas, arriving at a 
dista nce of 2.3 8 + " lTkpc. These distance measurements agree 
with lSato et al.l d2008l) within 2<x of the joint uncertainty. Also, 
our proper mot ion agrees within lcr with the proper motion by 

ISatoet al.ld2008l) . and our radial velocity (vlsr 29 km s ') 

is close to that of the water masers (vlsr ~ -31km s ). 
It is unlikely that the water and methanol masers originate 
at such different depths into the SFR that it would affect the 
distance. For example, in the SFR W3(OH), the water and 
methanol maser parallaxes are in good agreement dXu et al.l 
120061: lHachisuka et alJl2006l) . 

NGC281-W is an SFR th at lies near the edge of a su- 
per b ubble in the Perseus arm dMegeath et al.l 12003: Sat o et all 
2007). Our slightly shorter distance would change the position of 
NGC281-W on this super bubble. Another SFR located on this 
bubble has a par allax determined distance: IRAS 00420+5530 at 
2.17 ± 0.05 kpc dMoellenbrock et alj|2009l) . This allows to com- 
pare the distances and peculiar motions of both SFRs, and check 
whether the expanding super bubble, with the expansion center 
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Distance along 1=123 [kpc] 



NGC 281-W (Sato et al.) 



NGC 281-W (this work) 




IRAS 00420+5530 



- in 



12 3 

X [kpc] 

Fig. 21. Lower panel: a face-on view of the Galactic plane. 
The Galactic longitudes of NGC 281-W (/ - 123°) and IRAS 
00425530 (I - 122°) are indicated by dashed lines. The top panel 
shows a cross-section of the Galactic plane along the Galactic 
longitude of NGC 281-W (123°). The vertical coordinate of this 
plot is the direction perpendicular to the Galactic plane. We plot 
NGC 281-W and IRAS 00420+5530 as black dots, together with 
the respective distance error. Black arrows show the peculiar mo- 
tion of the source after removing the Galactic rotation (Rq = 8.4 
kpc, ©o = 254 km s ). The errors in peculiar motion are indi- 
cated by the gray cones. The light-gray circles mark the super 
bubble with the center of the bubb le at a distance of 2.5 kpc and 
a radius of 0.4 kpc, adopted after ISato et al.l d2008). The cross 
marks the averaged distance of 2.58 kpc. The dotted lines show 
the backward prolongation of the proper motion of NGC 281 - W 
and IRAS 00420+5530. 



at 2.5 kpc adopted after lSato et al.l d2008l) . can be responsible for 
the peculiar motions of both of them . We calculated the peculiar 
motions of IRAS 0042 +5530 after iMoellenbrock et al] 120091) 
and NGC 281-W after ISato etaT] d2008l) and our results. The 
bottom panel of Fig.l21lshows the peculiar motions in a face-on 
view of the Galactic plane. From the figure can be seen that the 
U component of NGC 281 - W, representing the vector toward the 
Galactic ce nter, is posi t ive an d non-zero in both our study and 
the work of ISato et all d2008l) . This means that NGC 281-W is 
moving toward the Galactic center, which is the expected move- 
ment for the near side of an expanding bubble. In the top panel of 
Fig.|2T|the peculiar motions of both SFRs are shown projected 
on a cross-section of the Galactic plane along the longitude 
of NGC 281 -W, I = 1 2 3°. Th e average distance to NGC 281- 
W based on ISato et all d2008l) and this work, 2.58+° g kpc, is 
marked in this panel as a cross. If we prolong the peculiar mo- 
tion backwards of both IRAS 00420+5530 and the average re- 
sult for NGC 281-W, they intersect at a distance from the Sun 
of 2.65 kpc and at a Galactic latitude of — 0°1. This suggets 
that, if indeed both SFRs are expanding linearly from one ex- 



pansion center, the expansion center is offset from the previously 
assumed position. However, this offset center of expansion is dif- 
ficult to understand when looking at the peculiar motion of both 
SFRs in the Galactic plane (Fig. [21] bottom panel), since the mo- 
tion does not seem to originate in one mutual expansion center; 
hence, it is difficult to pinpoint the characteristics of the super 
bubble if the peculiar motions and distances of only two SFRs 
are all that are known. 



5.6. S255 

For S255, the parallax distance, 1.59^}^ kpc, is much closer 
than the commonly used photom etric distance of 2.5 kpc 
dMoffat et al.lfl979l: iBlitz et al.|[l982h . S 255 is an individual Hn 
region, associated with a complex of Hn regions. The methanol 
maser emission coincides with a filament of cold dust and molec- 
ular g as between two Hn regions, S 255 and S 257. lMinier et al.l 
(2007) studied the star formation in this filament, which they 
propose is possibly triggered by the compression of the filament 
by the two Hn regions. They find several molecular clumps in 
this filament. At distance of m 2.5 kpc, the masses of the clumps, 
determin ed from the submi llimeter dust continuum, are around 
300 M dMinier et alj|2007l) . However, if we place S 255 at the 
parallax determined distance of 1 .6 kpc, the clump masses would 
drop by 60% to ~ 120 M . 



6. Summary 

We measured parallaxes of 6.7 GHz methanol masers using 
the EVN towards five SFRs, achieving accuracies approaching 
20 /vas. The primary results are summarized as follows. 

1 . We report trigonometric parallaxes for five star-forming re- 
gions, the distances to these sources are 2.57^ ^ kpc for 



ONI 



°' 776 -0 083 k P C fOT L 1206 ' °- 929 -0 033 k P C fOT L 1287 ' 



2.38+°;}|kpc for NGC 281-W, and 1.59^ kpc for S 255. 

2. Galactic star-forming regions lag circular rotation on average 
by ~17 km s _1 , a value comparable to those found recently 
by similar studies dReid et al.ll2009bl) . 

3. Individual 6.7 GHz methanol maser spots are stable over a 
period of > 2 years for most of the maser spots. The internal 
motions of the maser spots are weak (~ 0.5 - 1 km s ) and 
rectilinear. 

4. Measurements at 6.7 GHz are less disturbed by the tropo- 
sphere, as expected. However, the ionospheric delay cannot 
be ignored and is likely not to be completely removed by us- 
ing the JPL GPS-IONEX maps to calculate propagation de- 
lays. Continuum measurements of background sources show 
that most of them have significant structure that is evident 
from the large (up to 1-2 mas, or 2-10 km s -1 ) apparent 
movements between pairs of background sources. This ad- 
ditional uncertainty from the apparent movement does not 
prohibit a determination of the peculiar motion; however, it 
increases the error bars of the peculiar motion by the order of 
the uncertainty, depending on the distance and longitude of 
the source. The parallax uncertainty can be affected as well, 
if the cause of the apparent movement between the back- 
ground sources are internal structure changes, which do not 
need to be perfectly linear. Two background sources should 
be a minimum for astrometric measurements at 6.7 GHz, so 
three or more background sources are recommended. 
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Table 5. Detailed results of parallax and proper motions measurements 



Background 


V'LSR 


Parallax 


Ha 


Hs 


source 


(km s~') 


(mas) 


(mas yr~') 


(mas yr -1 ) 


ONI 






Northern group 




J2003+3034 


-0.4 


0.521 ±0.055 


-3.33 ±0.18 


-5.26 ± 0.41 




A A 


0.401 ± 0.052 


-3.72 ±0.17 


—5.26 ± 0.21 




0.4 


0.361 ± 0.052 


-3.82 ±0.17 


-5 28 + 19 




7 


0.526 ± 0.057 


-3.36 ±0.19 


_S i Q _i_ o 




Combined fit 


0.299 ±0.112' 








Avprapnnp* data 


0.391 ±0.061 








n a 
-u.*+ 


0.380± 0.121 


-2.87 ±0.39 


A QQ -i_ O ^1 
— H-.oa ± \J.J i 




u.u 


0.459 ±0.140 


-2.68 ± 0.46 


A OA . A 
-4.0U ± W.ZZ 




0.4 


0.483 ±0.142 


-2.59 ± 0.48 


-4.79 ± 0.20 




0.7 


0.415 ±0.112 


-2.81 ±0.37 


-4.78 ± 0.23 




Combined fit 


0.482 ±0.137' 








Averaging data 


0.368 ± 0.070 










Southern group 






1 A A 

14.4 




-3.88 ±0.12 


< m _t_ a 1 a 




14.8 




-3.75 ± 0.09 


-6.03 ±0.11 




15.1 




-3.67 ±0.21 


-6.09 ±0.19 


J2009+3034 


14.4 




-3.00 ±0.15 


-5.49 ± 0.20 




14.8 




-2.88 ±0.21 


-5.59 ±0.26 




15.1 




-2.77 ± 0.29 


-5.65 ± 0.35 


Both QSOs 


Combined fit 


0.390 ±0.1 16J. 








Averaging data 


0.389 ± 0.045 






< y" > north (o")pm 






-3.15±0.89(0.44£ 


-5.03 ±0.46(0.22) 2 


< H > south {CT)pm 






-3.33 ±0.90(0.45) 2 


-5.80 ±0.46(0.23) 2 


L1206 


J2223+6249 


-13.3 


1.163 ±0.222 


0.24 ± 0.85 


-2.84 ±0.40 




-12.9 


1.1 16 ±0.263 


0.08 ± 0.92 


-2.59 ± 0.48 




-10.9 


1.083 ±0.170 


0.23 ± 0.90 


-2.85 ± 0.23 




-10.5 


1.485 ±0.190 


0.56 ±0.31 


-1.61 ± 1.42 




Combined fit 


1.318 ±0.282' 








Averaging data 


1.331 ±0.180 






J2225+6411 


-13.3 


1.311 ±0.408 


0.22 ±0.51 


-0.28 ± 0.63 




-12.9 


1.322 ±0.386 


0.41 ± 0.50 


-0.65 ± 0.56 




-10.9 


1.300 ±0.416 


0.35 ± 0.52 


-0.50 ± 0.65 




-10.5 


1.174 ±0.237 


0.05 ± 0.64 


0.11 ±0.27 




Combined fit 


1.272 ±0.384' 








Averaging data 


1.288 ±0.241 






Both QSOs 


Combined fit 


1.331 ±0.250' 








Averaging data 


1.289 ±0.153 






< r 1 > (°")pm 






0.27 ±0.23(0.16) 2 


-1.40 ± 1.95(1. 15) 2 


L1287 


J0035+6130 


-27.0 


1.111 ±0.074 


-0.18 ±0.10 


-2.30 ±0.25 




-23.9 


0.928 ± 0.078 


-1.02 ±0.12 


-2.28 ± 0.08 




-23.5 


0.957 ± 0.093 


-0.93 ±0.11 


-2.48 ±0.12 




-23.2 


1.002 ±0.084 


-0.88 ±0.11 


-2.39 ± 0.09 




-22.8 


0.940 ± 0.083 


-1.14 ± 0.11 


-2.78 ±0.10 




-22.5 


0.917 ±0.067 


-1.17 ± 0.12 


-3.28 ±0.07 



1 The error of the combined fit multiplied by V5v, where N is the number of maser spots. 

2 We calculated an unweighted arithmetic mean of the individual proper motion results from all maser spots and background sources. The error 
bar on the mean is the standard error of the mean to which was added, in quadrature, the apparent movement between the two background sources 
of the respective coordinate. The uncertainty in the proper motion, which is introduced by the background source, is hereby taken into account. In 
parenthesis is given the standard deviation of the mean. 
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Table 5. continued. 



Background 


VLSR 


Parallax 


Ma 


Us 


source 


(km s" 1 ) 


(mas) 


(mas yr _1 ) 


(mas yr~') 




I nrYininpn Tit 
v^UlllUlllCLl 11L 


QRA + 086 1 










1 016 + OS? 






J0037+6236 


-27.0 


1.306 ± 0.07 1 


-0.14 ± 0.10 


-1.71 ± 0.40 




-23.9 


l.244 ± 0.044 


-0.88 ± 0.04 


-1.72 ±0.27 




-23.5 


1.225 + 0.040 


-0.81 ±0.04 


-1.77 ±0.25 




-23.2 


l.Oll ± o.ioi 


-0.88 ±0.14 


-1.83 ±0.12 




-22.8 


0.945 ±0.125 


-1.17 ±0.20 


-2.20 ±0.14 




99 S 




1 01 -+- O Ozt 


9 (\A. -+- O XI 

— Z.D+ ± V.J 1 




Combined fit 


1.192 ±0.107' 








Averaging data 


1.150 ±0.052 






Both QSOs 


Combined fit 


1.079 ±0.069' 








Averaging data 


1.077 ± 0.039 






< j" > (o-) P m 






-0.86 ± 0.1 1(0.33) 2 


-2.29 ± 0.56(0.46) 2 


NGC281-W 


J0047+5657 


-30.2 


0.416 ± 0.054 


-2.64 ±0.15 


-1.75 ±0.10 




-29.9 


0.380 ± 0.045 


-2.54 ±0.12 


-1.72 ±0.08 




-29.5 


0.401 ± 0.041 


-2.58 ±0.19 


-1.72 ±0.07 




-29.2 


0.404 ± 0.047 


-2.69 ± 0.22 


-1.75 ±0.08 




-28.8 


0.497 ± 0.026 


-2.55 ± 0.04 


-1.69 ±0.13 




-28.1 


0.529 ± 0.028 


-2.68 ± 0.04 


-1.52 ±0.14 




Combined fit 


0.400 ± 0.070' 








Averaging data 


0.398 ± 0.042 






J0052+5703 


-30.2 


0.459 ± 0.049 


-2.79 ± 0.06 


-1.87 ±0.27 




-29.9 


0.420 ± 0.048 


-2.68 ± 0.05 


-1.84 ±0.25 




-29.5 


0.388 ± 0.048 


-2.75 ± 0.05 


-1.84 ±0.22 




-29.2 


0.408 ± 0.054' 


-2.87 ± 0.07 


-1.87 ±0.20 




-28.8 


0.602 ± 0.056 


-2.65 ± 0.07 


-1.82 ±0.26 




-28.1 


0.648 ± 0.055 


-2.76 ± 0.07 


-1.65 ±0.34 




Combined fit 


0.399 ± 0.054' 








Averaging data 


0.425 ± 0.024 






Both QSOs 


Combined fit 


0.412 ±0.045' 








Averaging data 


0.421 ± 0.022 






<H> (cr) pm 






-2.69±0.16(0.10) 2 


-1.77 ±0.11(0.10) 2 


S255 


J0613+1708 


4.6 


0.628 ± 0.027 


-0.14 ±0.05 


-0.84 ± 1.67 








-0.14±0.54 2 


-0.84 ± 1.76 2 



